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ABSTRACT 

Hyper-accreting discs occur in compact-object mergers and in collapsed cores of mas- 
sive stars. They power the central engine of 7-ray bursts in most scenarios. We cal- 
culate the microphysical dissipation (the viscosity and resistivity) of plasma in these 
discs, and discuss the implications for their global structure and evolution. At the tem- 
peratures (kbT > m e c 2 ) and densities (p ~ 10 9 — 10 12 gr cm" 3 ) characteristic of the 
neutrino-cooled innermost regions, the viscosity is provided mainly by mildly degener- 
ate electrons, while the resistivity is modified from the Spitzer value due to the effects 
of both relativity and degeneracy. Under these conditions the magnetic Reynolds num- 
ber is very large (R c m ~ 10 19 ) and the plasma behaves as an almost ideal magneto- 
hydrodynamic (MHD) fluid. Among the possible non-ideal MHD effects the Hall term 
is relatively the most important, while the magnetic Prandtl number, Pm (the ratio of 
viscosity to resistivity), is typically larger than unity: 10 ^S Pm £6x 10 3 . Inspection 
of the outer radiatively inefficient regions indicates similar properties, with magnetic 
Prandtl numbers as high as ~ 10 4 . Numerical simulations of the magneto-rotational 
instability (MRI) indicate that the saturation level and angular momentum trans- 
port efficiency may be greatly enhanced at high Prandtl numbers. If this behaviour 
persists in the presence of a strong Hall effect we would expect that hyper-accreting 
discs should be strongly magnetised and highly variable. The expulsion of magnetic 
field that cannot be dissipated at small scales may also favour a magnetic outflow. 
We note that there are limited similaries between hyper-accreting discs and X-ray 
binary discs - which also have a high magnetic Prandtl number close to the black 
hole - which suggests that a comparison between late-time activity in 7-ray bursts 
and X-ray binary accretion states may be fruitful. More generally, our results imply 
that the possibly different character of high Prandtl number MHD flows needs to be 
considered in studies and numerical simulations of hyper-accreting discs. 
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1 INTRODUCTION 

The fluid dynamics of dense accretion discs is described 
by the Navier-Stokes and magneto-hydrodynamic equations. 
Within such a description, the microscopic properties of 
the fluid are parameterized (in part) via the kinematic vis- 
cosity v, which describes the rate at which inter-particle 
collisions damp fluid motions, and the resistivity 77, which 
quantifies the dissipation of currents via Ohmic losses. 
Timescales can be associated with each of these quanti- 
ties. The viscous timescale of a disc t acc = r 2 /u, for ex- 
ample, is the timescale on which the microscopic viscosity 



would lead to angul ar momentum transport over a scale 
r. As is well known (jPringlei Il981n . for astrophysical discs 
the timescales for evolution driven by the microscopic vis- 
cosity are many orders of magnitude larger than the ob- 
served evolution timescales. It is therefore generally assumed 
that evolution is driven by an "anomalous" or "turbulent" 
viscosity V\ ^> v l|Shakura fe Sunvaevlll973r ). possibly aug- 
mented in some cases by a ngular momentum loss in a wind 
JBlandford fe Pavng|l982l ). In most astrophysical discs v t 
is likely generated from turbulen ce driven by the m a gneto - 
rotational instability (MRI; see iBalbus fc Hawlevl (|l998h 
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and references therein). Since i>x S> v most studies of discs 
make the implicit assumption that the actual values of v and 
r\ are unimportant, and consider instead evolution under the 
action of the effective transport coefficients v t and r\t- 

Provided that the Reynolds (-Re) and magnetic 
Reynolds (R c m) numbers are large enough, the absolute 
value of either one may indeed be of little import for ac- 
cretion discs. Their ratio, however, which we define as the 
magnetic Prandtl number, 



Pm 



(1) 



is almost certainly not ignorable (jBalbus fc Hawlevl ll998). 
An accretion disc serves to convert the large-scale energy 
present in the shear flow into other energy forms - tur- 
bulent kinetic energy, magnetic field, and heat - and the 
dissipation required to generate heat may vary strongly 
depending upon the microphysical transport coefficients. 
In particular, if Pm S> I the plasma is highly viscous 
on the scale at which resistivity operates. In this regime 
one may expect that dissipation of the magnetic field will 
be suppressed, with a resulting modification in the sat- 
uration level of dynamo generated fields a nd an inverse 
casca de of magnetic energy to large scales IjBrandenburgi 
I2001T ). The importance of the magnetic Prandtl num- 
ber for the resulting magnetic field structure has been 
illustarted analytically ( Umurhan . Menou fc Regevl l2007al ; 



Umurhan, Rcgcv & Menou 2007b) and numerically, both I 
idealised simulations of dynamo a ction within forced tur- 
bulence (|Schekochihin et al.l 12004 ) and in local shearing- 
box simulations of the MRI within accr etion discs 
(|Lesur fc Longaretti|[2007l ; iFromang et al.ll2007l fl. For discs 
the current limited set of simulations indicates that the 
efficiency of angular momentum transport, parameterized 
via the Shakura-Sunyaev a parameter, increases with Pm, 
and can approach a ~ 1 given the combination of mod- 
est Pm = 8 and a w eak net vertical magne tic field. Moti- 
vated by these results. iBalbus fc Henril (2008'l calculated the 
expected radial variation of Pm within standard Shakura- 
Sunyaev models of geometrically thin, radiatively efficient 
accretion discs. They found that Pm exceeds unity within 
~ 50 Schwarzschild radii of the central object, and suggested 
that the qualitatively different behaviour of the MRI in the 
high Pm limit might be responsible for the time-dependent 
outbursts and state changes that are observed in X-ray bi- 
naries. 

In this paper, we compute the magnitude and ra- 
dial dependence of the magnetic Prandtl number in hyper- 
accreting discs. Such discs, which can support accretion 
rates of the order of 1 Mq s -1 , for m from compact object 
mergers (e.g. iRuffert fc Jankalll999T ) a nd when the co res of 
rapidly rotating massive stars collapse i|Wooslevll 19931 ). The 
disc structure is shaped by neutrino emission in the regions 
close to the accreting object, while in the outer parts neither 



1 Since dissipation in discs typically occurs on small scales for 
which the fluid is unaware of the large scale shear, one would 
anticipate that the basic physics of large Pm dynamos is similar. 
Disc simulations are required, however, to assess the impact of 
the small-scale physics on MRI saturation and the large scale 
flow properties. It is unclear whether the shearing-box approach 
is adequate for study of these effects. 



neutrinos (the temperature is too low) nor photons (which 
are trapped within the flow) can provide cooling. Hyper- 
accreting discs are invoked to provide the initial energy in- 
jection in many models for 7-ray bursts (GRBs). Our ulti- 
mate motivation for studying the magnetic Prandtl number 
in these discs is to understand the complex phenomenology 
(jets, highly variable prompt emission, X-ray flares, plateau 
phas es etc) that is observed in GRBs (e.g. iBurrows et al.l 
I2005T ) and which may well derive from the evolution of the 
disc. 

We anticip at e diff erent results from those found by 
IBalbus fc Henril (|2008f ) for photon-cooled discs. In the 
neutrino-cooled region the differences arise from two rea- 
sons. First, since neutrino opacities are vastly smaller than 
photon opacities, hyper-accreting discs are much cooler 
and denser than an extrapolation of Shakura-Sunyaev discs 
would imply, and this on its own would alter the magnetic 
Prandtl number. Second, the temperatures and densities in 
hyper-accreting discs are such that the electrons are mildly 
relativistic and mildly degenerate, while the elastic colli- 
sion cross-section for the nuclei exceeds the Coulomb c ross- 
section. As a consequence, the classical ISpitzerl |l962j) ex- 
pressions for v and 77, which are adequate for ordinary ac- 
cretion flows, no longer apply. They do apply, instead, in the 
outer regions of hyper-accreting discs, where a different Pm 
behaviour occurs since the flow is radiatively inefficient and 
photon-pressure dominated. 

The organization of this paper is as follows. In §2 
through §5 we describe the formalism for calculating the 
structure of hyper-accreting discs. This is required here in 
order to determine consistently the plasma conditions within 
these discs - including the temperature, density, nuclear 
composition and degree of degeneracy. Readers familiar with 
vertically averaged hyper-accreting disc models will find that 
our approach generally follows standard practice. The struc- 
ture of the resulting disc models is presented in §6. Given 
these conditions, we then calculate in §7 and §8 the resis- 
itivity and viscosity in the neutrino cooled inner region of 
the disc. These Sections contain the principle new results of 
this paper. Combining the viscosity and resistivity, we show 
in §10 that Pm > 1 or > 1 across the entire radial range 
relevant for GRB models. In §11 we estimate the strength 
of other non- ideal MHD effects, and in §12 we study the 
structure of the outer non-radiative region of the disc. §13 
and §14 discuss and summarize our results, and what they 
may imply for disc evolution and GRB observables. 



2 THE DISC MODEL 

We calculate the temperature and density at the midplane 
of a hyper-accreting disc from vertically averaged thin 
disc solutions. In hyper-accreting discs the diffusion time 
for photons to leak out from the disc is much larger than 
the accretion time scale (tdiff/tacc ~ 10 11 ) and photons 
remain trapped. However, for sufficiently high accretion 
rate^J the temperature in the innermost region exceeds 
~ 1 MeV and neutrino production switches on (Popham, 



2 IChen fc Beloborodovl 1120071 ) estimate a threshold value of 
M ;i 0.07 Mq sec -1 for a non-rotating black hole and a = 0.1. 
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Woosley & Fryer 1999). This happens for radii r ,$ 100 rs, 
where rs = 2GMbh/c 2 is the Schwarzschild radius for a 
black hole of mass Mbh. The neutrino flux liberates the 
energy deposited locally by viscous (turbulent) stresses and 
consequently the local sound speed c s remains smaller than 
the rotational velocity (v,/,) of the disc. Such discs are geo- 
metrically thin (the vertical scale height is small compared 
to the radius), and can be described by a modified version of 
the standard "a— disc" theory (Shakura & Sunyaev 1973). 
Many previous studies of hyper-accreting discs have a d opted 



this formalism (e.g. Popham. Woosley fc Frverl 



Narayan, Piran & Kumar 2001 ; Kohri fc M incshigc 



1999 



2002 



Pi Matteo. Perna fc Naravanl |2002| ; IChen fc Beloborodovl 
20071 . and references therein). In what follows we consider 
only the innermost region (r ^ 70 rs) where neutrino 
cooling is generally efficient. For simplicity the accretion 
rate is assumed to be time-independent, but our basic 
results for the magnetic Prandtl number would carry over 
locally to real discs in which the accretion rate varies with 
radius and time. 

In a steady state, conservation of mass and angular mo- 
mentum for a differentially rotating fluid yield a relation 
between the surface density E and the radial mass flow M, 



Eu t 



M_ 
37 



m 



1/2 



(2) 



where the fluid equations have been averaged in the ver- 
tical (z) direction over the pressure scale height H — 
—P/(dP/dz). We assume that no viscous torque is acting 
at the location of the innermost stable orbit r* = 3 rs. Con- 
sistent with the thin disc approximation radial pressure gra- 
dients have been neglected and the orbital velocity is simply 
i>0 = GMbh /r (the "Keplerian velocity"). The surface den- 
sity is related to the midplane (z = 0) density via E = pH. 
The radial mass flow M [gr s _1 ] is defined by the continuity 
equation 



M 



-27rrEt; r = constant, 



(3) 



where the radial drift velocity v r is sub sonic. The turbulent 
kinematic viscosity is parameterize d as (|Shakura fc Sunvaevl 
ll973l ; lLvnden-Bell fc Pringielll974r ) , 



v x = ac s H. 



(4) 



We emphasise here that the kinematic viscosity vx is quite 
different in nature from the "microscopic" (also referred to 
as "molecular") viscosity u, that enters into the definition of 
the magnetic Prandtl number. The first is associated with 
turbulent flow and magnetic stresses in the plasma, while 
the second is given by deflection of ions and electrons by 
particle-particle collisions in the plasma. 

For a thin disk both radial advection of energy and flux 
in the radial direction can typically be neglected (we justify 
this for our case later). From energy conservation the heating 
dissipated locally by viscous stresses is emitted per unit area 
through the surface of the disc as neutrino F v and anti- 
neutrino flux Fv : 



1 1/ E D. 2 = F v 

8 



F v . 



(5) 



The calculation of the flux in the optically thin and optically 
thick regions of the disc is detailed in the next section. Eq.[2] 
and eq. [5] allow us to solve for the disc structure once the 



equation of state is specified. In the inner regions that we are 
considering, densities (p ~ 10 8 — 10 11 gr cm -3 ) are high for 
the temperatures of interest (T ~ a few 10 10 K). This has a 
number of consequences. First, electrons are mildly degener- 
ate (j,/ 6 ~ a few, where p, is the chemical potential in units 
of the electron rest mass energy rn c c 2 and 8 = k^T/(m e c 2 ), 
where fcb is the Boltzmann constant. Second, the plasma is 
neutron rich. If Y p is the number fraction of protons over 
baryons, Y p ~ 0.1. Finally, the baryon pressure 



Pb 



k b T, 



(6) 



dominates the total pressure P (e.g. lBeloborodovll2003h . and 
the isothermal sound speed c 2 — P/p is, 



(m * 1/2 



cm/s. 



(7) 



At larger radii, on the other hand, the density is smaller, 
p/0 ^ 1, an d electron and photon pressure are more im- 
portant (e.g. iNaravan. Piran fc Kumarll200ll ). Note that in 
eq.[6] and henceforth, we ignore the mass difference between 
protons and neutrons. We also assume that all the baryons 
are free, since a particles present at larger radii have been 
photo-dissociated. This is a valid approximation since, for 
the parameters we use in our disc structure models, a par- 
ticl es do not constitute more t han 10% of the baryons (see 
e.g. IChen fc Beloborodovll2007l . thereafter CB07). 

The neglect of advection in eq. [5] is generally a good 
approximation. The ratio of neutrino diffusion time to ac- 
cretion time is 



L-acc C r 

~ 6 x 10 



"WJ I 0.3 



-/rs 
10 



-1/2 



(8) 



and hence advection is important only if the optical depth 
to neutrinos is of the order of 10 3 or higher. In our models 
optical depths are, at most, t„ ~ a few tens, unless the 
accretion rate is very large, (e.g. M Si 10 Mq sec -1 , for a — 
0.1). For a maximally rotating Kerr black hole (not treated 
here), advection can become important for lower accretion 
rates: M > 2 M sec _1 , for a - 0.1 (CB07). Advection of 
energy in photons instead can play a role in the region we 
consider for a ^, 0.1, since the outer radiatively inefficient 
part of the disc can extend inward, within 70 rs. We will 
comment on this in § [6] 



3 THE NEUTRINO FLUX 

At the temperatures and densities of interest for the inner 
regions of hyper-accreting discs neutrinos are produced pri- 
marily by electron capture onto protons 



p + e 



(9) 



while anti-neutrinos are produced by positron capture onto 
neutrons 



n + 6 1 



■ p + v. 



(10) 



The inverse of these processes - namely neutrino absorption 
by nuclei - constitute one source of neutrino opacity, others 
(discussed in § [5| are neutrino scattering by both nuclei and 
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electrons. We note that at threshold reaction [9] can proceed 
if the electron has an energy equal to the mass difference 
between a neutron and a proton. In this case the emitted 
neutrino has zero energjo The mild degeneracy of electrons 
suppresses the pair annihilation channel e + + e~ — » v + v, 
since positrons are exponentially reduced with respect to 
electrons: 



oc exp 



2// 



< 1. 



(11) 



At equilibrium and neglecting advection, the capture rates 
[51 and [TU1 should be equal (see also § Q. We therefore set 
F v — Fa and derive only the neutrino flux. The calculation 
basically follows the CB07 prescription. 

For e~ capture onto protons (O, the emissivity in units 
of m c c 2 is, 



e v = kp p \ (e - QYF(e, fi- )e\ / 1 - — de 



(12) 



(e.g. Shapiro & Teukolsky 1983), where e is the electron 
energy in units of rn c c 2 , Q — (m n — m p )/m e , the matter 
density of protons is defined as p p = Y p p, and k — 6.5 x 
10 -4 is a constant. The function F(e,p) is the Fermi-Dirac 
distribution, 



F(e,n) = 



exp[(^ii)] + r 



(13) 



The energy density in units of m e c 2 for an optically thin 
disc can be expressed in terms of the emissivity as 

U v , n = e„ — . (14) 

c 

If, conversely, the optical depth t„ 3> 1, neutrinos thermalize 
to a Fermi-Dirac distribution feq. H3|l and the energy density 
in units of m c c 2 is given by, 



U v ,c— — -V / F{e v ,n v )e v de v . 

27T 2 ft 



(15) 



For fx v — 0, U VtC — 24.5 k h T i . To deal with regions of the 
disc that have t v ~ 1 we interpolate between the optically 
thin and thick limits. We define 



C^,n 



Uu,n T U v ,< 



(16) 



and express the neutrino flux for arbitrary optical thickness 

as 



m c c 2 



U v , n c(l + T„) 1 for x < |, 

i7„,cc(l + r„)- 1 fora;>i 



(17) 



4 NUCLEAR COMPOSITION 

In the mid-plane of the disc pairs, baryons (free protons 
and neutrons only in our case) and photons are in thermo- 
dynamic equilibrium. For given (p, T) , and specified optical 
depth, this fact suffices to determine the nuclear composi- 
tion of the plasma. In practice we compute the composition 

3 A positron, instead, can be captured by a neutron at zero en- 
ergy, yielding a neutrino with non-zero energy. In the non degener- 
ate limit, this asymmetry favours an equilibrium shifted towards 
proton richness (see §|4]|. 



in the optically thick and optically thin limits separately, 
and deal with the intermediate regime by making an inter- 
polation that is consistent with charge conservation. 

Let us first note some general considerations. In addi- 
tion to the processes of electron / positron capture onto 
nuclei (equations l9l and I10|) . pairs are created and absorbed 



e + e ^7 + 7. 
This implies that 



(18) 



(19) 



Baryons, instead, are created and destroyed only according 
to the processes given by equations [9] and 1101 Therefore, 
at equilibrium (dY p /dt — 0), and neglecting advection, the 
emissivity for neutrinos and anti-neutrinos must everywhere 
be equal. 

In regions of the disc that are very optically thick to 
neutrinos, they will have relaxed into a Fermi-Dirac distri- 
bution and the above condition can be expressed as, 

n v = n 9 , (20) 

where the number densities are calculated by integrating 
eg. 1 131 over momenta with the statistical weight set to unity. 
The equality of the neutrino and anti-neutrino number den- 
sities (eq. I20[) implies a vanishing neutrino chemical poten- 
tial, 



Pv = p,v — 0. 



(21) 



The chemical balance of reactions [9] and \W\ (and their in- 
verses) fip + fi = fj, n + fi„ gives 



fi^ein— + Q, 



(22) 



l|Beloborodovll2003f ) where protons and neutrons (with num- 
ber density n p , n n respectively) have a Maxwellian distribu- 
tion. 

For regions of the disk that are optically thin, we use 

n v —hv-, (23) 

where the number of neutrinos v and anti-neutrinos v emit- 
ted per unit volume per unit time is 



n v = kp p \ (e - Q) F(e,p-)e^\jl ^ de, 



and 



(24) 



ho—kpn (e + Q) 2 F(e,p+)e 2 Jl -de, (25) 

(e.g. Shapiro & Teukolsky 1983). 

In both regimes, we impose the conditions of charge 
neutrality, 



n- + n+ = Y p , 

m p 

where 
v — n P 

1 V 3 

n b 
and baryon conservation, 
n„ 

1 D- 



Tip ~\~ Tin 



(26) 

(27) 
(28) 
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The electron n_ and positron n+ number densities are ob- 
tained by directly integrating eg. 1 131 over momenta with the 
statistical weight set to two. To summarise, we use two dif- 
ferent sets of equations in the two extreme regimes: when 
t u ^> 1, we find Y p = Y PjC and p, — p c from eqs. 1221 [26l 
and 1281 when r„ <C 1, we find Y p = Y Pin and p = p n from 
eqs. [231 [261 and [281 

Matter in the disc plane, of course, can have an arbi- 
trary optically thickness. To compute the nuclear composi- 
tion for specified p and T we first interpolate to find the 
chemical potential using, 



jJb— (1 — x)p n + Xp c , 



(29) 



where x is given by eq. 1161 We then find the proton fraction 
Y p by imposing charge conservation (eq. [26}. The resulting 



Jp IS ^p,n ^ Ip 5^ Jp,c- 



Physically, the most important trends in the nuclear 
composition can be understood in terms of the ratio of the 
temperature to the (relativistic) degeneracy temperature 



Tdc 



0.71 



pYp 



1/3 



K, 



(30) 



where we set n p = n_ = pYp/m p . If T < Td cg (at low 
temperature, or at high density) then p, increases, positrons 
are suppressed feo.lll[). and the resulting plasma is neutron- 
rich (low Yp). In the opposite limit of T > Tdcg, n+ ~ 
n 
energetically favoured. 



but Yp H 0.5 because, as noted before, reaction [10] is 



5 OPACITIES 

The vertical optical thickness of the disk to v and v in eq. 1171 

r v = Hpn, (31) 

is calculated considering three sources of opacity: k = K a + 
k s , c + K Sj b [cm 2 gr -1 ]. First, absorption onto nuclei, 

«a = — K, n (l - Y p ) + a^pYp) (32) 



10 - -. 







10 



0.2 M sec 

o 

2 M sec" 1 



4 M sec -1 



10 



r/r„ 



30 



60 



Figure 1. Midplane temperature as a function of radius (in units 
of the Schwarzschild radius rg), for three disc models computed 
for Mbh = 3Mq and a = 0.1. The accretion rates are indicated 
in the legend. 



Then, the usual interpolation allows us to link the optically 
thin and thick regimes: 



E v — (1 - x) Ev,n + X -E„,c, 

(CB07) where 

P - £v 
n v 

corresponding to r„ C 1 and 

k b T 



E 



3.17 



(35) 



(36) 



(37) 



corresponding to t„ ^> 1 is the mean energy calculated from 
a Fermi-Dirac distribution function with zero chemical po- 
tential. We use the same treatment for anti-neutrino cross 
sections. 

In practice, we find that for conditions appropriate to 
our disc models absorption by nuclei is always the domi- 
nant opacity by one or two orders of magnitude. Electron 
scattering is least important. 



where o"„. n ( p ) is the cross section for v (ii) absorption by 
neutrons (protons). Second, elastic scattering by electrons, 



k SiC = 0.5— (a c>u + er ^) . 



(33) 



where Y c — n^/n^ and cr eiV tp\ is the electron-neutrino (anti- 
neutrino) cross section. We neglect scattering onto positrons 
since n+/n_ <C 1. Finally, elastic scattering by nuclei, 

k, = (0.3 Fp + 0.36 (1-Y p ))^-El, (34) 

m p 

where ctq = 1.7 x 10~ 44 cm 2 . For the explicit forms of the 
cross sections, we refer the reader to the Appendix B of 
CB07. 

The cross sections for these processes depend on the 
neutrino energy, and the total cross section a should be de- 
rived by integrating over the neutrino spectrum. This re- 
quires knowledge of the neutrino distribution function at 
arbitrary optical depth. We, instead, use approximate mean 
cross sections, which are adequate for our purposes. We first 
evaluate a at a mean neutrino energy (in unit of m e c 2 ) : 

a = (j{E v ). 



6 RESULTS FOR THE DISC STRUCTURE 

To give an idea of the midplane temperatures and densities 
encountered in hyper-accreting discs, we show in Figure \T\ 
and Figure [2] illustrative examples of the radial profiles of 
these quantities. The disc models plotted were computed 
for a black hole of 3 Mq, with a = 0.1 and accretion rates 
of 0.2, 2 and 4 MQsec -1 . For these parameters the disc is 
optically thin across the entire radial range shown for the 
lowest value of the accretion rate. For the two high accretion 
rate models the disc is optically thick within about 15rs of 
the black hole. 

The most important feature of the neutrino cooled so- 
lutions for our purposes is the relatively weak radial depen- 
dence of the temperature. As noted elsewhere (e.g. CB07), 
the disc self-adjusts to maintain electrons in a state of 
mild degeneracy (i.e. T w 0.3 — 0.5Td cg ): higher degener- 
acy would cause the cooling to decrease drastically and the 
temperature to rise, lifting the degeneracy, and vice versa 
for lower degeneracy. This thermostatic aspect of neutrino 
cooled discs has as a consequence that temperature does 
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Figure 2. The same as fig. [T] but for the midplane total matter 
density. 



not vary strongly with radius, since the degeneracy temper- 
ature depends only weakly on density (and in addition Y p 
decreases with p). The proton fraction Y p is everywhere less 
than 0.3, falling to around Y p ~ 0.1 in the innermost regions. 
We have verified that the physical conditions in our disc 
models reproduce those of prior authors. Comparing our re- 
sults to those of CB07, we find excellent agreement - espe- 
cially for a < 0.1 - between our models with a given a and 
their model with a/2 (the factor two arises primarily due to 
different definitions of viscosity and different inner boundary 
conditions, both of which are largely arbitrary). For a = 0.1, 
their densities are lower and the advection dominated outer 
region (which is not cooled efficiently by neutrinos) extends 
to smaller radii. Since we ignore the advective term in the 
energy equation, our profiles have a different slope but they 
never differ by more than a factor of 2 in density and by 
50% in temperature. 



V = 



7b 



4irr B n-t c 



(38) 



(see e.g. Jackson 1963 p. 460, for the non-relativistic deriva- 
tion). In this expression n_ is strictly equal to n-(p): the 
number density of electrons with momentum p, 7 C is the 
electron Lorentz factor, r c — e 2 /m c c 2 ~ 2.82 x 10 _13 cm is 
the classical electron radius and t cp is the proton-electron 
deflection time. The associated electrical conductivity is 
a — c /(t]4it). Extending the non-relativistic definition (that 
can be found for example in Longair 1992, p. 304), t cp can 
be defined so that /Apj_) t ep = p 2 , where (Apj) is the 
mean square impulse received by an electron per scattering 
perpendicular to its direction of motion. We obtain 



p 2 v 



87rn p e 4 logA p(p) ' 



(39) 



where logA cp (p) is the Coulomb logarithm, defined as the 
logarithm of the maximum to the minimum impact param- 
eter for proton-electron collisions. 

If we take the non-relativistic limit (^ e — > 1) of eqs. [35] 
and [39] and insert the mean thermal v elocity of non- 
degenerate electrons v = v c — y/ 'Sk-^T '/ 'ivm c , we recover the 
scalingsj derived by Spitzer (1962) for a completely ionized 
plasma with n p = n_ 



r, s = 1.9xl0 12 logA ep , s T- 3 / 2 ^L 

sec 



(40) 



where the Coulomb logarithm for electron-proton collision 



A 



o P ,S 



1.24 x 10 4 



T 3/2 1 4.2 x 10 5 



(41) 



In this non-degenerate regime, the (inverse) dependency of 
7/s on density is weak. 



7 ELECTRICAL RESISTIVITY OF THE 
PLASMA 

The computation of the resistivity of a plasma in which elec- 
trons are both relativistic and degenerate is quite subtle, 
and ultimately requires numerical methods. We first review 
some general considerations, including the classical Spitzer 
(1962) result, which we emphasise is not applicable to hyper- 
accreting discs. These results will prove useful later when we 
discuss the relative importance of different non-ideal MHD 
effects. 



7.1 General considerations 

The electrical resistivity of a plasma is a measure of the 
electron mobility in the presence of an applied electric field. 
Acceleration due to the electric field (eE, where e is now 
the electric charge) is opposed by drag from (primarily) 
electron-proton collisions. Balancing these competing forces 
on a electron of velocity v and momentum p = y e m e v, 



eE - 



u 



where we neglect pressure gradients and magnetic fields and 
assume stationary protons. Ohm's law then gives, 



7.2 Numerical computation of r\ 

To compute the plasma resistivity under our condit i ons o f 
tem perature and density, w e use the code bv lPotekhinl (|1999h 
and lPotekhin et al.l (119991 1 5 I. which is able to treat arbitrary 
degeneracy of fully relativistic electrons (but not positrons) . 
The calculation involves an integration of eq. [38] over en- 
ergy, weighted by the Fermi-Dirac distribution feq. H3[) . The 
Coulomb logarithm is calculated via an effective scattering 
potential that takes into account Debye charge screening 
(see reference above for details). We use as input param- 
eters to the code the temperature T and the proton den- 
sity p p , which, in turn, depends (via Y p ) on T and p. From 
p p the code derives the electron density and thus the elec- 
tron chemical potential which uniquely specifies the electron 
Fermi-Dirac distribution. 

In principle the calculation of the total resistivity should 
include the contribution of positrons both as charge carriers 
(they contribute to the current density and have their own 
resistivity) and as deflectors for electrons. However, we have 
already noted that positrons are exponentially suppressed 
with respect to electrons and from charge conservation it 



4 There is a numerical factor of 3.5 between the Spitzer resistivity 
and our derivation. 

5 http://www.ioffe.ru/astro/conduct/ 
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follows that they can also be neglected in the drag force 
which remains dominated by protons. 

The code allow us to explore the plasma resistivity un- 
der different condition of electron degeneracy. In the non- 
degenerate limit the resistivity decreases with temperature, 
because scattering becomes less frequent, while the density 
dependence is weak and confined to the Coulomb logarithm 
(e.g. Spitzer resistivity, eq. 1401 for k^T <C m c c 2 ). In the 
highly-degenerate regime, instead, the resistivity increases 
with temperature, because more energy states become avail- 
able for the electrons to be scattered into. The density de- 
pendence here is strong. Higher density (for fixed tempera- 
ture) results in higher degeneracy and suppression of scatter- 
ing: the plasma resistivity thus decreases. The intermediate 
mildly degenerate regime, where the resistivity changes very 
slowly with temperature and density, is the one occupied by 
our plasma. 



7.3 Results 

The calculated resistivity is shown as a function of tem- 
perature in Figure [3] For a density of p = 10 10 gr cm~ 3 



the degeneracy temperature Td cg: io (evaluated at p p = 
0.5 x 10 10 gr cm -3 ) is approximately 10 11 K. For tempera- 
tures above this value, the proton fraction remains approx- 
imately constant at a value of Y p ~ 0.5, thus p p does not 
change appreciably and the resistivity decreases with tem- 
perature since electron-proton collisions becomes less fre- 
quent. For lower temperatures, one's naive expectation is 
that the resistivity ought to flatten out and eventually de- 
crease as electrons become more and more degenerate. This 
behaviour is seen if we artificially impose a temperature in- 
dependent electron density (dashed line), but it does not re- 
flect the physical behaviour of the disc plasma. In fact, as the 
temperature drops toward and below Td eg ,io the plasma be- 
comes more and more neutron-rich (i.e. Y p decreases), with 
an attendant decrease in the electron density. The effective 
degeneracy temperature decreases, and the electrons never 
become highly degenerate. The overall result is that at "low" 
temperatures the electrons self-adjust toward a mildly de- 
generate regime in which 77 rises slowly with decreasing tem- 
perature. The dependence on density is weak. 

Within a hyper-accreting disc, the midplane tempera- 
ture is a fairly weak function of radius (Figure[T|. This, when 
coupled with the weak dependence of the resistivity on den- 
sity, results in a radial profile of r/, shown in Figure [4] that 
increases slowly with distance from the black hole. Between 
5 and 70 rs the variation in the resistivity is less than an or- 
der of magnitude. There is a modest decrease with increasing 
accretion rate. 



8 MICROSCOPIC VISCOSITY OF THE 
PLASMA 

Contributions to the viscosity of the plasma come from three 
sources, proton-proton Coulomb collisions, elastic baryonic 
scattering due to the strong force, and electron viscosity. It 
is not immediately obvious which of these processes dom- 
inates for a mildly relativistic and degenerate neutron-rich 
plasma. We first consider plasma viscosity due to baryon 
collisions. Typically (but not in our case) only ion-ion 



10 
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10" 



10 
Temperature [K] 

Figure 3. The microscopic viscosity (upper curves) and electrical 
resistivity (lower curves) as a function of temperature, for three 
different densities: p = 10 9 gr cm -3 (dot-dashed line), p = 10 10 
gr cm -3 (solid line) and p = 10 11 gr cm -3 (dashed line). These 
curves are calculated with a self-consistent nuclear composition 
and chemical potential p. The resistivity curve plotted with a 
dashed line with dots is computed for p = 10 10 but with a fixed 
nuclear composition, Yp ~ 0.5, and the corresponding p for each 
temperature: this forces the electrons to become degenerate as T 
decreases. The proton fraction Y p ~ 0.5 corresponds to the maxi- 
mum temperature (~2x 10 11 K) considered here, for any shown 
density. All curves are computed in the optically thin limit. The 
vertical lines show the degeneracy temperature (eq. 1301 . calcu- 
lated for a given density with Y p ~ 0.5. 
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Figure 4. Microscopic viscosity (upper curves) and electrical re- 
sistivity (lower curves) as a function of scaled radius, for the three 
disc models shown in Figure [T] 



Coulomb collisions need to be taken into account. Protons 
are non-degenerate and non-relativistic, so we use the non- 
relativistic limit of eq. I39I with p = m p v where the mean 
proton relative velocity v = w rB i = V2v is writ ten in terms 
of the mean thermal velocity v = ^JSk^T/ivm-p. Multiplying 
by (log A eP: s/log Ap P ,s) we obtain the proton-proton (p-p) 
deflection timescale, 



ipp.s = 18.03 



T 



3/2 



n p logA P p, s 



(42) 



where the p-p Coulomb logarithm is given by 
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10 

Temperature [K] 

Figure 5. Viscosity components as a function of temperature for 
matter density p = 10 10 gr cm -3 . The corresponding curves are 
labelled. The electron viscosity v e appropriate for the partially 
degenerate regime is shown as the light dashed line. This compo- 
nent, which is the one that we use in the disc calculation, domi- 
nates the total viscosity (shown as the solid curve) at low temper- 
atures. At high temperatures the viscosity due to non-degenerate 
electrons (t'e.ndegi shown as the heavy dashed line on the right) 
increases, but in this regime it is always negligible compared to 
the contributions from neutron- neutron (u n ) and proton- neutron 
(i>p) collisions. At the upper left of the plot we show the value 
of the electron viscosity in the strongly degenerate regime (which 
we emphasize is not appropriate for disc conditions). 



A pp ,s = 1.24 x 10' 



rp3/2 



(43) 



(|Spitzerl 1 1962T l 6 l. Scattering in this regime is dominated by 
small-angle deflections. 

Once plasma temperatures exceed ~ MeV, however, 
as in hyper-accreting discs, the elastic nuclear cross sec- 
tion becomes larger than the Coulomb one. The shape of 
the nuclear two-body cross section can be derived theoreti- 
cally with a partial-wave analysis (considering only s- waves) 
withi n the frame-work of the effective-range theory (JBethei 
1 19491 ). For the particle energies of interest to us, the par- 
tial cross section (for a single total spin state) depends only 
on two parameters: the scattering length a and the effective 
range of the nuclear potential r ef , 



°0,i = 



47ra; 



[1 - |a ; r c f,iK 2 ] 2 + (oik) 2 



(44) 



where the wave number is k — (l/h)y^2m p E, where E — 
(3/4)fct>T is the centre-of-mass total kinetic energy. For the 
p-n collision we need to consider contributions from both 
singlet i — s (antiparallel spins) and triplet i — t (parallel 
spins) states, and the total cross section is the weighted sum 
of partial siglet and triplet cross sections, 



= -ffo.t + 7 Co,. 



(45) 



with at = 5.41 x 10~ 13 cm, a s - -23.71 x 10~ " cm, r ef . t = 
1.74 x 1CT 13 cm and r el>s = 2.75 x 1CT 13 cm (|Hackenbursl 



2006). For identical particles, however, the triplet state is 
excluded. The p-p and n-n cross section are thus, 

1 



'pp 



tCTq, 



S! 



(46) 



with a s = -17.31 x 10"- 
Cr nn = (l/4)<70.s, with a, 



cm and r c f >s 

= -18.8 x 10~ 13 cm and r rf|S = 
2.75 x 10~ 13 cm (|Miller et alj|l990h . The result is that a pn > 

Cnn — Cpp • 

Once we allow for the possibility of strong interactions 
the total collision timescale for protons becomes 



1 



+ ■ 



+ 



1 



(47) 



where the nuclear p-p and p-n timescales are t PPinu = 
l/(°p P npi>rei) and t P n = l/(cr pn n n v TC \) respectively. The col- 
lision timescale for a neutron is instead 

I = JL + -L K JL, (48) 

t t i t 

t-n <-n P fcnn tiin 

where t np = l/(a pn n p v s ) and t nn = l/(cr nn n n v s ). Note that 
for both protons and neutrons the collision timescale is dom- 
inated by encounters with neutrons, since the flow is neutron 
rich and the high temperature inhibits electromagnetic col- 
lisions. The corresponding mean free paths are A p = vt p and 
A n = vt n . 

Beside baryons, we should consider the contribution to 
viscosity from electrons, since degeneracy results in a longer 
mean free path A c = (3 ct ap (J3 is the ratio of electron velocity 
to c) and a higher momentum density than non-degenerate 
electrons. 

According to the classical kinetic theory the dynamical 
viscosity r\ v [gr cm -1 sec -1 ] for a single species "i" is propor- 
tional to pi^iAi. Therefore, we can write 

t]v = C\ {p P vX p + pnVXn + n-p\ c ) (49) 

where C\ is a constant of the order of unitjlj. The typical 
electron kinetic energy is fim c c 2 , where /i is the c hemical 
potential, thus its momentum is p = 
The total kinematic viscosity v 
written as 



;v V + 2/i. 
T)vl p [cm 2 s x \ can be 



;^s 



v ^P 
P X~s 



Y, 



A n t n-p A c 
pv A s 



p; A^ + 



(50) 



where pus — CipvXs an d As = vt pp ,s (where t pp ,s is calcu- 
lated for n p — n_ = n) is the dynamical viscosity for a fluid 
of completely ionized hydrogen, 



u s = 2.21 x 10" 



rpb/2 



plogA pp ,s 



(51) 



(Braginskii 1957, 1958; Spitzer 1962). In eq. [50] the first 
term on the right is the proton viscosity u p , the second 
is the neutron viscosity u p and the third is the electron 
viscosity u c . This last reduces formally to the formula 
for strongly degenerate elect rons pf c ,dc g oc n— pf vf t cp 
l|Nandkumar fc Pethickl|l984| ) (with Ci = 1/5). 

In fig. [SJ we plot v and its components as a fucntion 
of temperature. In the non-degenerate limit (T > Tdeg), 



Our definition differs by a factor of 
Spitzer result. 



1.6 from the classic 



7 We assume that the constant C\ is the same for the collision 
processes considered here. 
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protons and neutrons contribute equally to the total vis- 
cosity: they carry the same momentum density, transported 
in both cases by n-p collisions. In this regime, the elec- 
tron contribution f e ,ndeg has been calculated as !/ e ,ndeg = 
vs \Y e v% t Ep /v \ C p,s) , where Y B is the electron fraction and 
t cp = 0.15 sec T 3 ' 2 /n p log A eP: s (eq. [39] with v — v e and 
p = v c m c ). As the temperature decreases, the flow be- 
comes neutron rich and the electrons degenerate. As this 
happens the proton momentum density drops faster than 
the neutron's one, strong collisions with neutrons dominate 
the baryon collision frequency and the electron mean free 
path increases (while the baryon one keeps decreasing). As 
a consequence, electron viscosity becomes comparable or 
greater than the baryon viscosity. Since the electron density 
decreases with temperature, v c eventually decreases (even 
if A c keeps increasing) and it does not attain the limiting 
viscosity value for strongly degenerate electrons (shown in 
the figure as the horizontal dashed line). Therefore the total 
viscosity v decreases as the plasma cools even when electron 
viscosity dominates. The behaviour of v as a function of 
density is shown in Figure ED It decreases as the plasma be- 
comes denser, since p-n collisions become more frequent and 
the electron (and thus their momentum) density decreases 
when below the nominal degeneracy temperature. In these 
figures, the T < Tdcg regime represents the typical situation 
in hyper-accreting discs. 

Figure [4] shows the predicted radial profile of the viscos- 
ity within the hyper-accreting disc models. In the optically 
thin regions of the disc the viscosity increases with radius, 
since the temperature in this region decreases much more 
slowly with radius than the density. Electron viscosity is 
higher than baryon viscosity but the two contributions tend 
to converge (and they can become comparable for low ac- 
cretion rates) in the less dense regions away from the black 
hole where the electrons are less degenerate. The trend is 
the opposite in the optically thick innermost region, which is 
present only at high accretion rates. In this case the viscosity 
decreases moving outward. Here, the neutron contribution 
is highly suppressed (~ one-two orders of magnitude) and 
electrons completely dominate the flow viscosity. 



8.1 Neutrino and photon viscosities 

Photons and neutrinos can also in principle contribute to 
viscosity. In the case of neutrinos, the mean free path is of 
course enormous compared to that of the baryons, so mo- 
mentum transport by neutrinos will dominate if it acts as a 
viscosity. This will depend upon the neutrino optical depth. 
In regions of the disc that are very optically thick neutrino 
viscosit y will act as a (very large) Navier-Stokes viscosity. 
Indeed, Masada et alj (|2007 ) have suggested that neutrinos 
might be able to shut off the MRI completely, by reducing 
the Reynolds number below that which admits linear growth 
of MRI mode^j. In optically thin regions, on the other hand, 
momentum transport by neutrinos is no n-local and acts as 
a radiation drag rather than a viscosity l|Arav fc Beeelmanl 
1 1992 ). In this paper we are principally concerned either with 



8 Note that one cannot easily appeal to neutrinos to shut off 
turbulence entirely, since the neutrino flux itself arises as a con- 
sequence of heating driven by the very same turbulence. 
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Figure 6. Comparison between the "X-ray binary" disc model 
and the "neutrino cooled" disc model for the electrical resistivity 
and the microscopic viscosity as a function of scaled radius. The 
"X-ray binary" disc has a IOMq central object, M = 0.1L e dd/c 2 
(where the Eddington luminosity is ~ 10 39 erg s _1 ) and a = 0.01. 
The neutrino cooled disc model has M = 0.2MQsec — 1 , a = 0.1 
and M BH = 3 M Q . 



optically thin regions of the disc (where neutrinos do not 
contribute viscosity), or with regions of modest optical depth 
[t v ~ a few) for which the correct treatment is unclear. The 
reader should, however, bear in mind that the presence of 
neutrinos may have dynamical consequences in the optically 
thick regions of the flow. Since we derive high values of the 
Prandtl number while ignoring the neutrino contribution, to 
the extent that neutrinos contribute to viscosity, they will 
only reinforce our conclusions on the large Pm regime of 
hyper-accreting discs. Our magnetic Prandtl numbers are 
conservative in that sense. 

The discussion on the different regimes according to the 
optical thickness also applies to photon viscosity. Our discs 
are everywhere in the mid-plane extremely optically thick 
to photons. This implies a much smaller mean free path for 
photons than for neutrinos and their contribution to viscos- 
ity is relatively suppressed by a factor equal to the ratio of 
mean free paths for the two species. Photon viscosity can 
also be neglected with respect to baryon viscosity, even if 
they have similar mean free paths, since the photon mo- 
mentum density is small compared to the baryons' one. 



9 COMPARISON WITH X-RAY BINARY 
DISCS 

Identical methods can be used to compute the resistivity and 
viscosity for photon-cooled discs, which are present in X-ray 
binaries at similar radii around black holes of a few Solar 
masses. In Figure [S] we plot the radial profiles of these quan- 
tities for a model presented by Balbus & Henri (2007) as rep- 
resentative of a X-ray binary disc. The model has a central 
black hole of mass 10 Mq, surrounded by a disc accreting at 
one tenth the Eddington limit with a — 0.01. We treat the 
photon flux in the diffusion approximation, adopt an opacity 
that is of Kramers' form plus an electron scattering correc- 
tion (k ~ 0.43 + 6.6 pT~ 7 ' 2 ), and write the total pressure 
as the sum of the gas pressure (baryons and electrons) and 
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Figure 7. Magnetic Prandtl number as a function of scaled radius 
computed for the three neutrino cooled disc models shown in 
Figure [T] The magnetic Prandtl number for the "X-ray binary" 
disc model is also shown for comparison, as a dotted curve. 




Figure 8. Contours of magnetic Prandtl number in hyper- 
accreting discs as a function of radius and accretion rate, com- 
puted for a = 0.01. 



photon pressure. The abundances are taken as 10% a par- 
ticles and 90% hydrogen by number. The gas is completely 
ionised. In the X-ray binary regime the Spitzer resistivity 
(eq. I40p is valid (we neglect e — a collisions) , while the vis- 
cosity is given by p-p Coulo mb collisions and a-p C oulomb 
collisions (see Appendix A in lBalbus fc Henri 2008). 

In the X-ray binary disc model the plasma is cooler and 
less dense, and electrons are non-degenerate. These charac- 
teristics enhance the e-p collision probability, resulting in a 
resistivity that is several orders of magnitude higher than in 
neutrino cooled discs. The viscosity is also much higher, and 
moreover, it decreases towards larger radii. This is the op- 
posite trend to that seen for the viscosity in neutrino cooled 
discs. 



10 THE MAGNETIC PRANDTL NUMBER 

With the results for the resistivity and viscosity in hand we 
can compute the magnetic Prandtl number, defined as the 
ratio of the microscopic viscosity to the resistivity, 

Pm = -. (52) 

n 

Results for the three disc models that we have discussed in 
detail are shown in Figure [7] The radial behaviour of Pm is 
mainly determined by the viscosity, since the resistivity is a 
rather weak function of radius. In optically thin regions of 
the disc (everywhere for M = 0.2 MQsec -1 , and at large 
radii, r ^ 15 — 30 rs, for the higher accretion rate models) 
the Pm number is relatively flat but eventually increases 
with radius because the rapidly decreasing density increases 
the microscopic viscosity in the disc. In the optically thick 
regions, the temperature increases more rapidly toward the 
black hole than does the density, and the increased viscos- 
ity leads to larger Pm in the central parts. We note that, 
although the values of both the resistivity and the viscosity 
are far larger in the photon cooled X-ray binary disc model, 
the value of the magnetic Prandtl number is not all that 
different in these inner parts. More significant is the differ- 
ent radial behaviour. For hyper-accreting discs Pm increases 
with radius once the disc becomes optically thin, and this 




Figure 9. The same as Figure. [8] but for a = 0.1. 



trend continues throughout the neutrino cooled zone. For 
X-ray binary discs Pm falls off with radius and eventually 
becomes smaller than one (at r ~ 45 r g ). O ur results for this 
model are qualitatively similar to those of iBalbus fc Henril 
(2008) (see their fig. 1, left panel) with the quantitative dif- 
ferences arising due to different boundary conditions and 
value of r]s ■ 

Figures [8] [9] and \W\ show the dependence of Pm on 
radius and accretion rate for three different values of a, 
a = 0.01,0.1 and 1. We find that Pm > 1 throughout the 
neutrino cooled region for all accretion rates in the range 
0.2 M e sec -1 < M < 2 M Q sec" 1 . Typically Pm > 1 (ex- 
cept in the r > 20rs regions for a = 0.01), with values as 
high as a few 10 3 if a — 1. 

As a decreases, the disc density is higher and the viscos- 
ity drops: the lower viscosity results in a smaller magnetic 
Prandtl number. This is unlike X-ray binary discs, where the 
a dependence of the microscopic vi scosity is largely offset by 
a similar change in the resistivity (|Balbus fc Henrill2008r ). 

All the results discussed here are computed assuming 
midplane disc conditions. Moving toward the upper layers 
of the disc the density invariably drops more rapidly than 
the temperature, and as a consequence the electrons become 
non-degenerate. However, verticle mixing is efficient and the 
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Figure 10. The same as Figure. [8] but for a = 1. 

nuclear compos ition (Y v ) at higher latitude is close t o that at 
the midplane (|Rossi. Armitage fc Pi Matteol 120071 ). Under 
these conditions the resistivity increases more slowly than 
the viscosity (whose main contribution will be eventually 
come from neutrons) and Pm increases further. 



11 OTHER NON-IDEAL MHD EFFECTS 

Our motivation for computing the magnetic Prandtl num- 
ber derives from various suggestions that the saturation level 
of the MRI varies strongly with Pm. Almost all numerical 
simulations addressing this issue include ohmic dissipation 
as the only non-ideal term in the MHD equations. If other 
non-ideal terms are dominant in hyper-accreting discs there 
will be (at the very least) additional uncertainties in inter- 
preting the results of simulations for astrophysical disao 
Here, we estimate the magnitude of the different non-ideal 
terms for fiducial values of temperature, T — 3 x 10 10 K, 



and density, p — 5 x 10 10 gr cm 3 



the plasma composition is Y p 



Under these conditions 
5.5 x 10 -2 , the chemical 



potential is fi — 8.7 and the sound speed c s = 1.6 x 10 cm 
s _1 . The equipartition magnetic field, 

B cq = V8^K (53) 

is 1.8 x 10 15 Gauss. 

The low resistivity of plasma in neutrino cooled discs 
results in a very large magnetic Reynolds number. If we 
take H = 5.3 x 10 6 cm (~ 0.3 x 20 r s ), then 

cJJ_ 
V 



RbM — 



1.5 x 10 1 



(54) 



where we have used the consistent value of the resistivity 
7/ = 5.85 x 10" 4 cm 2 s" 1 . All non-ideal MHD effects are 
therefore small. In relative terms, however, the high density 



9 Most simulations to date have examined the Pm 2> 1 regime 
in the absence of the Hall effect, while simulations that have in- 
cluded the Hall term have focused on the regime where the mag- 
netic Reynolds number is low (as is appropriate for protoplane- 
tary discs). Although the Hall term should not directly alter the 
dissipative properties of the fluid, we cannot rule out that the 
behavior of a high Prandl number disc with a strong Hall effect 
might differ from the pure Pm S> 1 case. 



means that the plasma is sufficiently collisional that, for 
very weak magnetic fields, ohmic dissipation is stronger than 
both the Hall effect (the separation of charge carriers in the 
presence of a magnetic field perpendicular to the current) 
and ambipolar diffusion (the differential motion of neutrons 
and charged particles) . The dominance of ohmic dissipation 
persists until the magnetic field becomes strong enough that 
the cyclotron frequencies of the charged particles exceed the 
collision frequencies. We can write the relat ive strength of 
the o hmic (O) and Hall (Ha) effects as (e.g. ISano fe Stone I 

EIxl), 

Ha , . 

— =— = oj cc £ C p. (55) 

where uj ca is the electron gyro-frequency and t cp the e-p 
deflection timescale. For our conditions 

fi + 1 



7e ^ 

47rr e ?7n_ 3.54 x 10~ r2 r]np 



2.87 x 10" ls sec, (56) 



where the electron density has been evaluated as ri- 
pYp/rrip, and 

1.76 x 10 7 B TT 
■ Hz, 



(7e> 



gauss 



(57) 



where the mean Lorentz factor is 7 C ~ /i + 1. The field B 
here is the mean local field. Eq. [5S]then gives, numerically, 



^=5.2X10--^. 
O gauss 



(58) 



Likewise, the ratio of ambipolar (A) to Hall (H) strength 
can be written as fe.g. ISano fc Stone 1 12002). 



A_ 
Ha 



(2 - 2Y P ) = 9.4 x 10" 



B 



gauss 



(59) 



is the proton gyro- frequency 



where lu cp = o; cc (7 e m c /r 
and the p-n collision time-scale is t p 
We note again that the Coulomb coupling (t cp ) is weaker 
than the strong collision coupling for our plasma conditions. 
Eq. [55] and eq. [59] show that ohmic dissipation is the 
dominant non-ideal MHD term only in weakly magnetised 
plasmas which have a ratio of thermal to magnetic energy 
P > 10 8 . This means that, for any reasonable magnetic field 
in the disc, the Hall effect is the dominant non-ideal MHD 
effect. For our fiducial parameters the Hall effect becomes 
the most important non-ideal MHD effect for B £ 2 x 10 11 
C1 1Q I. At equipartition ambipolar diffusion is also more im- 
portant than resistive diffusion but we remain in the Hall 
regime. 



12 THE RADIATIVELY INEFFICIENT ZONE 

In the outer disc — beyond ~ one hundred rs — the phys- 
ical conditions are quite different. In this region the disc is 
photon-pressure dominated. The temperature (T £, 10 10 K) 
is below the neutrino cooling threshold and the flow is ef- 
fectively non- radiative. The free baryons are in the form 
of H e nuclei IjChen fc Beloborodovl 12007! ) and the density 



10 When u> cc t cp Js 1, t cp and r\ are formally the relaxation 
timescale and the resistivity along the magnetic field lines. The 
transverse currents have a collision timescale reduced by 1/(1 + 
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(p ^ 10 7 grcm~ 3 ) is low enough that the electrons are non- 
degenerate (and non relativistic) . The disc is thus both more 
resistive and more viscous, and the resistivity and viscosity 
have the classical Spitzer dependencies on temperature and 
density (eq. [40] and eq. 15 1 p but with numerical coefficients 
(also in the "Spitzer" Coulomb logarithms eqs. 1411 and I43f) 
adjusted to be appropriate for a completely ionised He flow. 
It is worth emphasizing that unlike the radiatively iricffi- 
cien t flows discussed in the X-ray binary or AGN context s 
fe.g lMahadevan fc Quataertlll997l ; lTanaka fc MenoullJOOr}) . 
those in hyper-accreting scenarios would remain fully colli- 
sional. This constitutes a rare situation in which the colli- 
sionality guarantees a unique temperature for electrons and 
nuclei. Moreover MHD dissipation scale arguments based on 
classical Spitzer values are adequate. 

The magnetic Prandtl number in this region is, 



rp4t 

Pm = 7.3 x 10" 28 — -, 



M' 



(60) 



where A = logA C H C x logAHoHo is the product of the 
Coulomb logarithm for e — H c and for H c — H c collisions. 

If we assume that the temperature is set by a fraction 
of the virial energy density then aT 4 = / [GMhh/r)p. Nu- 
merically we find that T 4 /p = 5.9 x 10 34 / (r s /r) and 



Pm 



1.44 x 10 b 



(?) 



(A/15) I,/' <61) 

where we take / = 1/2 1_|. A varies very slowly with radius 
and eq. |6l]with (A/15) = 1 is actually a very good approxi- 
mation for Pm in the whole outer disc. We note that Pm ~S> 1 
and it falls off slowly as r _1 : for r Si lOOrs, Pm ,$ 10 4 . 

This result indicates that in the outer regions of hyper- 
accreting discs Pm has a very different behaviour than in the 
innermost parts. Here Pm is rather insensitive to the actual 
accreting conditions M and a (as long as the outer disc re- 
mains photon-pressure dominated), but it is a stronger func- 
tion of temperature and it decreases with radius. However, 
both regions have in common Pm > 1: unlike X-ray binary 
and AGN discs, hyper-accreting discs seem not to have a 
Pm < 1 region. Beside a high Pm, this outer disc flow has 
other characteristics in common with the neutrino-cooled re- 
gion. The conductivity is high with R e M ~ 10 12 — 10 11 and 
the Hall effect dominates even for modest magnetic fields 
{fl < 10 6 — 10 5 ). These features are caused by the still high 
temperature that makes Coulomb coupling rather inefficient 
(though in the outer region it is still more efficient than nu- 
clear coupling). 



13 DISCUSSION 

In this section, we attempt to relate the plasma properties 
that we have computed to the global properties of discs. 
There are many steps required to connect the microscopic 
scales where resistivity and viscosity act to the global evo- 
lution of the disc, and simulations and analytic arguments 
provide only an imperfect guide. However, plausible guesses 



11 We are aware that there is no consensus on the exact value of 
/; however, the dependence of Pm on / is small. 



are possible. We focus on the possible consequences of mag- 
netic Prandtl number greater than unity for the magnetisa- 
tion and evolution of hyper-accreting discs. 

Let us first summarise the argument that suggests that 
the magnetic Prandtl number matters. In turbulent plasmas 
with Pm > 1, the absence of turbulence below the viscous 
scale A„ ~ Re H means that the resistive scale may not be 
easily accessible to field lines and reconnection could be slow 
with respect to magnetic amplification. We can understand 
the dependence on Pm with the following simple argument. 
Turbulence at scale A„ exponentially amplifies the magnetic 
field on an eddy turnover time-scale £ e ddy = H/c s (\ V /H) ' , 
while the magnetic field dissipation on that same scale takes 
t v = \t/rj. The ratio of these timescales is 



*-cddy 



Pm. 



(62) 



So, for Pm > 1, magnetic energy can build up, first on the 
viscous scale and subsequently in an upward cascade toward 
larger scales l|Brandenburgil200H ). 

This may have important consequences for accretion 
discs. Within discs, the turbulent growth of magnetic fields 
can be limited - in principle - by either resistive dissi- 
pation (on small scales) or buoyant expulsion of the flux 
from the disc. Numerical simulations suggest that in low 
Pm discs resis tive dissipation ult imately sets the satura- 
tion level (e.g. IStone et al. II1996T ). If, however, the small- 
scale dissipation of magnetic fields is frustrated by a high 
magnetic Prandtl number, the only channel available to 
limit the flux may be macroscopic buoyancy. Since buoy- 
ancy is relatively inefficient - the timescale is likely to mod- 
estly exceed the Alfven crossing time H/Va of the disc - 
a high magnetic Prandtl number disc could saturate at a 
level that involves stronger magnetic fields, more vigorous 
turbulence and angular momentum transport, and ongo- 
ing field expulsion. These expe ctations are suppor t ed by 
the available numeri c al evi dence. iLesur fc Longarettil (|20070 
and iFromang et al.l 1)20071 ) find that the parameter a in- 
creases with the magnetic Prandtl number in the range 
2 < Pm < 8 an d R e ~ 10 3 - IP 4 . I n the case where a 
approaches unity, ILesur fc Longarettil (|2007i ) find a highly 
time dependent disc, with large fluctuations in a and the 
mean magnetic pressure. This appears to be a consequence 
of the well-known dependence of the strength of the MRI 
on the magnetic field: as the ratio j3 of thermal to mag- 
netic pressure decreases a at fir st increases until the MRI is 
quenched around equipartition l|Hawlev. Gammie fc Balbua 
ll995l ; lLesur fc LongarettHl2007l ). 

Based on these arguments, two evolutionary conse- 
quences can be contemplated for hyper-accreting discs 
formed after stellar core collapse or compact object merg- 
ers. In the initial phase, a high magnetic Prandtl number in 
the neutrino cooled inner disc results in stronger magnetic 
field and a higher efficiency of angular momentum transport. 
Since, as we have shown, a higher a in such discs results in 
a yet higher Pm, a runaway occurs until a ~ 1. This regime 
is sustainable for many dynamical timescales, provided that 
the excess magnetic energy that could quench the MRI is 
removed from the disc, for example, by buoyancy. This is 
plausible since close to equipartition the buoyancy timescale 
~ H/Va ~ O -1 is similar to the orbital timescale and can 
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provide cooling by removing magnetic energM 12 !. The disc, 
which would be in a sim ilar regime to that simulated by 
iLesur fc Longarettil (|2007l ). would be highly variable. This 
initial stage, in which a highly variable magnetic "corona" 
is built above the disc via flux expulsion, could be associated 
with the production of a magnetic outflow, d yna mically sim- 
ilar to those stu died by Proga et al.l (120031 ) and iMcKinnevI 
IJ200rJ ) (see also iKing et all 120041 ) , polluted with some disk 
matter, that eventually will be observed as a GRB. 

The longer term evolution of very high Pm discs is less 
clear. One possibility is that eventually a net vertical field, 
which cannot be removed locally by buoyancy, builds up and 
suppresses the MRI. This depends upon the ability of field 
lines to diffuse outwards. Turbulent diffusion can depend in 
turn on the strength of the magnetic field and on non-linear 
MHD effects. In our case the dominance of the Hall effect 
may influence the behaviour of the turbulence, by introduc- 
ing an asymmetry between magnetic fields that are aligned 
or anti- aligned with the dire ction of the angular momentum 
vector (|Sano fc Stone 1120021 ). The outcome is thus uncertain 
and the investigation requires numerical methods. 

From a purely phenomenological perspective, one might 
note that there are similarities between the state transitions 
in X-ray binaries and GRB prompt and late-time activity. 
Our results show the presence of a high Pm number inner 
region in both discs. One may then speculate that the GRB 
prompt phase could involve similar physical conditions, in- 
sofar the magnetic field behaviour is concerned, to the tran- 
sition to the "low/hard state" in X-ray binaries, which is 
observationally asso ciated with the emergence of an outflow 
l|Fender et al.lll9 99). However, the theoretical analogy with 
X-ray binaries is imperfect, since it is unclear that what 
causes a change in state may be the same physical mecha- 
nism in both sources. For X-ray binary discs the temperature 
dependence of Pm may allow a local thermal runaway, and 
it has been suggested that this ulti mately causes global state 
transitions (jBalbus fc HenriluOOa) . The conditions for such 
a mechanism to work in hyper-accreting discs are present 
only in the outer radiatively inefficient regions, since in the 
neutrino-cooled regions Pm is primarily a function of den- 
sity. However, a classical thermal instability is not possible 
in radiatively inefficient flows. 



14 CONCLUSIONS 

In this paper, we calculated the magnetic Prandtl number 
for plasma conditions in hyper-accreting discs. This dimen- 
sionless number Pm = u/rj — A^/A 2 expresses the relative 
size of two critical turbulent flow scales: the scale on which 
velocity fluctuations are viscously damped \ v , and the scale 
A I; on which ohmic losses result in magnetic field dissipation. 
We find that in the inner neutrino-cooled regime: 



12 Note that we investigated a disc structure with a = 1 in which 
a fixed fraction of the cooling occurs non-radiatively via loss of 
magnetic energy. We found that Pm is only reduced by a factor of 
~ 2 with respect to a disc with a = 1 in which only neutrino cool- 
ing and gas pressure are considered. Hence reasonable amounts 
of non-radiative cooling do not avert a runaway to a high a state. 



(1) Electric resistivity involves relativistic, mildly degener- 
ate electrons. As a consequence, the resistivity is very low 
and weakly dependent on density and temperature. 

(2) The main source of viscosity is Coulomb collisions be- 
tween mildly degenerate electrons and protons. It decreases 
as the flow gets denser or cooler, since the electron fraction 
decreases. 

(3) The magnetic Prandtl number is always greater than 
unity, unless the angular momentum transport efficiency, 
parameterized via the Shakura-Sunyaev a, is very small 
(q < 0.01). Pm ranges between a few tens to a few 10 3 , and 
it typically increases with a. In the optically thin regions 
for neutrinos, the main dependence of Pm is on density (in- 
versely, via viscosity) since the disc temperature varies very 
slowly as a function of M and radius. This causes Pm to de- 
crease with increasing accretion rates. In the optically thick 
regions, the main dependence is on temperature, that varies 
more rapidly than density. As a consequence, Pm increases 
with M 

In the outer, radiatively inefficient region: 

(4) Resistivity and viscosity have the classical "Spitzer" de- 
pendencies on temperature and density, since electrons are 
non-relativistic and non-degenerate and viscosity is mainly 
given by Coulomb interactions between ionised helium parti- 
cles. The magnetic Prandtl number is thus a much stronger 
function of temperature Pm oc T 4 /p. In contrast, it does 
not depend on the overall magnitude of accretion (via M 
and a). The Prandtl number is always much larger than 
unity (Pm ~ 10 4 (100 rs)/r) and it is unlikely to fall below 
unity for the entire extent of the disc. 

We conclude by noting a number of implications for the 
global structure of hyper-accreting discs: 

(5) For all plausible values of the magnetic field strength in 
the disc, the Hall effect provides the largest non-ideal term 
in the MHD equations. These discs are perhaps a unique 
example of a high magnetic Reynolds number flow in which 
the Hall effect is the dominant non-ideal term. Since there 
are no simulations where these conditions are accounted for, 
their effect on the magnetic field evolution and MRI is un- 
clear. 

(6) The large values of the magnetic Prandtl number mean 
that the evolution of magnetic field within hyper-accreting 
discs is likely to be qualitatively different from discs with 
lower magnetic Prandtl numbers. The results of forced dy- 
namo and MRI simulations suggest that, in the high Pm 
regime, small scale field dissipation is suppressed and the 
saturation level o f the magnetic field is enhanced, per- 
haps dramatically (ILesur fc Longarettill2007l : lFromang et al.l 
|2007| ). Numerical simulations of hyper-accreting discs that 
do not account for microphysical dissipation may therefore 
have underestimated the magnetic field strength. 

(7) The immediate consequence of a high value of Pm for 
models of the central engines of GRBs is that disc formation 
is likely to be accompanied by vigorous expulsion of mag- 
netic flux. This favours models in which energy is liberated 
in the form of a magnetic jet or outflow. 
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